THE DYNAMICS OF OFF-CENTER REFLECTION 
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Abstract. Dynamical properties of a two-parameter circle map, called off-center re- 
flection, are studied. Certain symmetry breaking phenomena in the bifurcation process 
are illustrated and discussed. 



Background and Definition 

We study the dynamics of a two-parameter family of circle maps R r n : S 1 — > S 1 called 
the off-center reflection. When Q = n, this map is a one- dimensional analog of the general 
map raised in [|Y|, problem 21], which is geometrically a reflection in the circle. For other 
values of Q, it can be seen as a reflection with a deviation between the reflected and 
incident angles. Iterations of this map are not the natural sucessive reflections in the 
circle; nevertheless, this map is interesting for various reasons. The off-center reflection 
can also be seen as a perturbation (with small r > 0) of rotation by Q. In fact, it has 



an analytical form which extends the well-known Arnold circle map, lArlf . Its dynamics 



is related to the perturbation properties of Mathieu type differential equation, [|Ar2 . 
Furthermore, when the perturbation parameter r goes to 1, the map goes to another 
famous circle map, the doubling map. 



The off-center reflection is introduced in |ALj] by the following geometric description. Fix 



a point L inside the unit circle S 1 . For a point <p G S 1 , a ray is emitted from L to 0. 
This ray is "reflected" to hit S 1 again at a point, denote R r ^{(j)) in the future. This 
point is defined to be the image of under the map. It is quoted "reflected" because 
the "reflected" angle has a constant deviation ZL ^ from the incident angle l(4>). To 
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have the map geometrically well-defined, there is some restriction on Q. However, we 
will see later that it is analytically meaningful for other Q. In fact, it is sufficient to 
consider Q e (—71,71}. Furthermore, since the action has certain symmetries, it is no loss 
of generality to assume the point source L at (r, 0) with < r < 1. This is why the 
off-center reflection is given by the two parameters (r, Q). 




In ||ALj| , particular interest is placed on R r ^, i-e., when the reflected angle equals the 
incident angle. There, the link between dynamics and contact geometry of the map is 
studied. Here, we deal with the dynamical properties of the family with some focus 
on fi = 0, 71. It is because in these two cases, the off-center reflection map repects the 
symmetry of the circle and symmetric cycles may occur. We are particularly interested in 
when symmetric cycles of the map break into asymmetric ones. The properties presented 
in this article may be considered the first steps to understand the dynamics of the map. 
It is expected that deeper studies may bring forth more understanding to general circle 
maps. 

In §1, we will introduce the basic analytical properties of the map. In §2, the attracting 
orbits, especially the symmetric ones, of the map are investigated. Then, finally in §3, the 
bifurcation of the map is looked into. In particular, we give an explanation of how and 



OFF-CENTER DYNAMICS 
when the symmetric orbits go through a period preserving pitch-fork bifurcation. 

analysis in []Br| of similar behavior among certain cubic polynomials is borrowed. 



3 

The 



1. Analysis 

For our purpose, we consider S 1 C C ~ R 2 and it is covered by R under the exponential 
map x i— > e lx . Then the off-center reflection 



i ► RrM) : S 1 — > S 1 



has a unique continuous lifting, R r ,n : R — > R, which takes to fi. Since i? r; n(^ + 27r) = 
-Rr,n(^) + 27r, we may focus our attention on the interval (— tt, tv]. Let the incident angle 
be denoted by l(x) for x G (— 7T, tt]. Since i(a;) — > as a; — >• ±7r, it defines a continuous 
27r-periodic function on R, which is also denoted as l(x). In fact, it can be written in 
terms of the principal argument (with values between — tc and 7r) as 

i(x) = Arg(cosx — r + isinx) — x; 

and R r n(x) = x + Q — 2i(x). The first few derivatives of R are listed below. 



l'{x) 



ricosx — r) 



(cos x — r) 2 + sin 2 x ' 



1 - 2t'(x) = - 

2r(l — r 2 ) sinx 



1 — 4r cos x + 3r 2 



cosx — r) 2 + sin 2 x 



'cos x — r) 2 + sin 2 x 



2r(l -r 2 ) 


(1 + r 2 ) cosx — 2r(l + sin 2 x) 




(cos x — r) 2 + sin 2 x 


3 
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oo r k 

This incident angle has a series expression given in ||AL , l(x) = ^ — sin (for). Therefore, 

fc=i k 

the lift of Rr t a mapping to Q is given by 

R r ,n( x ) = x + Q — 2l(x) 

oo r k 

= x + Q — 2 ^ — sin(foc). 
fc=i ^ 

From this, we see that the Arnold circle map x i— > x + fi — £ sinx can be seen as a truncated 
version of the off-center reflection. The technique of our analysis in this paper may also 
be adapted to a similar study on the Arnold circle map. 

We will first establish a nice property for the dynamics of a map. The dynamics of the off- 
center reflection in the range 0<r<l/3is relatively simple. The following proposition 
is useful for studying the range r > 1/3. 

Proposition 1.1. For 1/3 < r < 1 and all Q, R r & is a map with negative Schwarzian 
derivative. 

Proof. The Schwarzian derivative of R r ^i is given by 

R'" Q {x) 3 (R" Q (x)\ 2 r(l-r 2 )H(r,cosx) 



where 



R' r ^(x) 2\R' rn (x)) (l-4rcosx + 3r 2 ) 2 (l-2rcosx + r 2 ) 2 
H(r, cos x) — (2 + 28r 2 + 6r 4 ) cos x 



13 + 19r 2 — (1 — r 2 ) cos 2x + 4r cos 3x 



I4r + 18r 3 ) + (2 + 40r 2 + 6r 4 ) cosx 



+ 2r(l — r 2 ) cos 2 x — 16r 2 cos 3 x 



Let y = cosx and consider H{r,y) for y G [—1, 1]. We have 
H(r, -1) = -2(1 + r) 3 (l + 3r) 



-2 - 12r - 24r 2 - 20r 3 - 6r 4 < 0; 



i/(r, 1) = 2(1 - r) 3 (l - 3r) < 0, for 1/3 < r < 1. 



OFF-CENTER DYNAMICS 5 

Moreover, its derivative wrt y satisfies 

d y H = (2 + 40r 2 + 6r 4 ) + 4r(l - r 2 )y - 48r 2 y 2 ; 

and for 1/3 < r < 1, 

d y H(r, -1) = 2(1 + r)(l - r 2 )(l - 3r) < 0, 
d y H(r, 1) = 2(1 - r)(l - r 2 )(l + 3r) > 0. 

We see that, inside (—1,1), the cubic polynomial y t— > H(r,y) can only have a single 
critical point and it is a minimum. Hence H(r,y) < max{if(r, — l),H(r, 1)} < 0. □ 

2. Attracting Orbits 

It is easy to see from our earlier formula that R' rn (x) > for < r < |, with equality 
only when r = | and cosx = 1. Thus, R r> Q is a homeomorphism for 0<r<l/3, so 
the dynamics is trivial. On the other hand, R r ^ is only a degree 1 map for r > 1/3. We 
would like to explore the dynamics of it in the coming sections. 

In the study of periodic orbits of R r> n, the information about the function i is often 
helpful. Since i is a 27r-periodic odd function, it is sufficient to know its properties in the 
interval [0,7r]. To be precise, i[x) > and is concave down for x G (0,7r). Its maximum 
value of 7r/2 — a r is attained at a r where a r is the angle satisfying < a r < ir/ '2 and 
coset r = r. A picture of its graph will be helpful to see its properties. 
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Furthermore, when r varies from to 1, l(x) varies from the constant zero function to 
a discontinuous linear function. These will be useful in calculations involving periodic 
orbits of R r ,n- For example, one may determine the birth of fixed point according to this 
knowledge of i. 

Proposition 2.1. For all r and Vt £ [tt — 2a r , tt + 2a r ], the map R t q has no fixed point 
and a saddle-node bifurcation occurs at \VL — ir\ = 2a r , that is, when r is the cosine of the 
deviation of between incident and reflected angles. 

Proof. To find fixed points of R r ,n, one tries to solve the equation 

Rr,n(x) = x + Q — 2i(x) = x mod 2tt, 

which can be reduced to 

— = dx) mod tt. 

Since t(— a r ) < t(x) < i(a r ), for Q e (— tt, tt], the equation has solution if and only if 

— TT Q TT 

— + a r <-<--a r . 

The fixed point of i? rj n(x) at the boundary parameter values Q = ±(tt — 2a r ) occurs at 
x = a r . It is easy to see that for < r < 1 and any f2, R! rCl (a r ) = 1 and R" n (a r ) ^ 0, so 
as (r, Q) crosses this boundary, a saddle node bifurcation occurs. □ 

1 + 2r 2 

Corollary 2.2. Let a r as above and b r e (0, a r ) such that cos6 r = — . The region 

that R r ,n has attracting fixed point is 

{(r,Q) : 2t{b r ) < \Q\ < 2i(a r ) } . 

In fact, the equation 2i(b r ) = \Q\ determines the generic values of r for the happening of 
period- doubling bifurcation. 
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Proof. If x G (— 7r, 7r] corresponds to an attracting fixed point of -R r ,Q, we have R r ^(x) = x 



mod 27T and 
to 



|1 — 2t'(rr)| < 1. From the expression of l'(x), this is equivalent 



cos a r = r < cos a; < 



1 + 2r^ 



3r 



Then, x G (— a r , — 6 r ) U (6 r ,a r ). Since t(x) is increasing on both (— a r , — b r ) and (b r ,a r ), 
Q/2 lies in the intervals defined by the image under t, namely, t(6 r ) < |fi/2| < t(a r ). To 
see the period doubling bifurcation, it is sufficient to show 



dr 



7^ 0, along the curves = ±2t(6 r ). 

x=±b r 

Here we need some calculations which will also be useful in the future. 



Lemma . We have 



dR rt n(x) —2 sin a; 

dr 1 — 2r cos a: + r 2 ' 

dR' n (x) Ar - 2(1 + r 2 ) cosx 



dr 

Furthermore, 

d(R 2 r n)'(x) dRUx) ~, 



(1 — 2r cosx + r 2 ) 2 



dr 



^ . R' rn (R r ^( x )) + R' r Q ( x ) 



dEL 



dr 



R r ,ni x ) 



dR r ,n(x) 
dr 



In particular, if R r ^(x) = x or R r n(x) = —x, one has 
d(R 2 r ,n)'(x) ~, dR! r ^x) 



dr 



dr 



1 + 



dR r>n (x) 
dr 



Proof of the lemma. The first two results only require simple calculus. The third one is 



a repeated application of the chain rule. Then, using that both R' r n and 
functions in x, the last result follows. 



dR' r 



r,n 



dr 



are even 



□ 



At Q — ±i(b r ), we have R r ,n(b r ) = b r and R' rn (b r ) = —1, therefore 



dr 



d(R 



x=±b r 



dr 



1 + 



dR 



x=b r 



dr 



x=b r 



Moreover, for r > 1/2, 



OR' 



dr 

9R r n 
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4r — 2(1 + r 2 ) cos6 r 



-6(1 - 2r 2 ) 



x=b r 



dr 



(1 — 2r cos& r + r 2 ) 2 r(l — r 2 ) 
—2 sin b r 



-2^4r 2 - 1 



x=b r 



1 — 2r cos 6 r + r 2 



7^-1, 



except at r - ,, 
at Q = ±2t(b r ). 



15+V241 



0.5119. Thus, generically, period-doubling bifurcation occurs 



□ 



The above information about fixed points is illustrated by the following picture, in which 
the grey areas are the region in (r, f2)-plane where fixed points exist. The darker area is 
where attracting fixed points occur. 




We would like to look into maps in the family that respect the symmetry of the circle. 
Let p : S 1 — > S 1 be the reflection across the real axis. Then i? rj n o p = p o R r fl if and only 
if fl = 0, 7r. This symmetry corresponds to the fact that R r ^(x) — Q is an odd function. 

For Q — 0, 7T, if n is the smallest integer that R™ n (<f)) = p(<ft) ^ <ft, one can easily show that 
<p belongs to a periodic orbit of prime period 2n. We call it a symmetric orbit of period 2n. 
In terms of Rr^n, this corresponds to R™ n (x) = —x ^ x mod 2ir. An orbit is asymmetric 
of period 2n if i? 2 ^(0) = <fi but R^((f>) ^ p{4>)- If an asymmetric orbit is formed by <p, 
then another asymmetric one, called the twin orbit, is formed by p(4>). The twin orbit 
may be itself when p(0) = (p. For the off-center reflection, a self-twin asymmetric orbit 
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must be the 2-cycle |e l0 ,e 17r |. There are numerous articles about symmetric periodic 
orbits in dynamical systems, especially on continuous types. An early one is [0. Their 
attention is on the return map of some reversible mechanic system such as the three-body 
system. 

Proposition 2.3. For any r > 1/3 and Q = 0,tt, the map R r Q has either no attracting 
orbit, or one symmetric attracting orbit, or two (counting multiplicity) asymmetric twin 
attracting orbits. 



Sketch of proof . We mainly use two properties of R r ^ to conclude this. First, each R r ^ 
has negative Schwarzian derivative, the technique of [CE] or [Br | is applicable. Therefore, 
once the map R r ^ has an attracting cycle, at least one critical point of it must be attracted 
to this attracting orbit. Since R r ,a, has only two critical points, there are at most two 
attracting orbits. If the attracting orbit is a reflection symmetric one, then this orbit 
attracts both critical points. If the orbit is an asymmetric one, its twin orbit attracts 
another critical point. Multiplicity occurs if there is a self-twin orbit. □ 



Besides the above general information about symmetric and asymmetric cycles, it is in- 
teresting to know two specific cases about it. The first is about self-twin asymmetric 
2-cycles. 

Proposition 2.4. A single asymmetric 2-cycle for R r Q occurs if and only ifVt = n. The 
orbit is { e i0 , e m } which is attracting when r < 1/ a/5 . 



Proof. The orbit for a single asymmetric 2-cycle can be obtained by solving the equation 



10 THOMAS AU 

The only simultaneous solution exists when Q = n and = e l0 , e ln . Then, whether the 
orbit is attracting can be easily decided by calculating R! (0)R! (n). □ 

Secondly, one expects that symmetric cycles occur naturally for Q = 0, ir. Do such cycles 
exist even if the map is not "symmetric" ? The following result provides a partial evidence 
for the answer. 

Proposition 2.5. The map R r ^ has a symmetric 2-cycle if and only ifQ = 0, ir. In both 
cases, the symmetric 2-cycle is unique. The cycle of R r0 is attracting while that of R r ^ 
is repelling. 

Proof. Let |e 1:c ,e~ 1:r | be a symmetric 2-cycle, i.e., R rt n(e ±lx ) = e Tlx . Equivalently, 

x + — - = l[x) mod 7T, 



and also, 



x — — = l(x) mod 7r. 



If we subtract the second equation from the first one, we obtained that Q = mod n 
and £1 E (— 7r, it]. Clearly, = 0, n. 

For = 0, the trivial solutions x = 0, ir corresponds to fixed points of i? r ,o- Further 
calculation on i shows that the equation x = i(x) has only a nontrivial solution for 
r > 1/2. In fact, let Ci E (0, 7r/2) such that cosci = then ±c x form a symmetric 



2-cycle. Since R' r0 (ci)R' r0 (— c\) 



l-3r 2 n ' 



, it follows that the cycle is attracting for 



- < r < —=. Moreover, since \i{x)\ < vr/2, one can only obtain trivial solution 0, it from 
2 y2 

x = i(x) + lor for k 7^ 0. 

For fl = 7r, the above equations always has solution for all r. In fact, its solution within 



1 _ _|_ 8r 2 

(— 7r, 7r) is given by the following. Let c 2 E (7r/2,7r) satisfy cosc 2 = . We 
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claim that C2 and 2tt — C2 are solutions to the equation x — | = i(x) and x + | = 
respectively. Taking tangent to both sides and using that i(x) = Arg(cos x — r + i sin x) — 
Arg(cosx + isinx), we have 

— cos x r sin x 
sin x 1 — r cos x 

It follows that 2r cos 2 x — cos x — r = 0. Hence j e 1C2 , e _1C2 j is a symmetric 2-cycle. It can 
be easily checked that 



Therefore, the symmetric 2-cycle is repelling. Again, by that \l(x)\ < 7r/2, we do not have 
other solutions to x ± | = + A;7r for A; 7^ 0. □ 



3. Bifurcation 



In this last section, our aim is to understand bifurcations between symmetric and asym- 
metric orbits of this family of maps. In particular, we would like to address some of the 
questions on the dynamics of R r ^ raised in [AL] . Let us first look at the asymptotic orbit 
diagrams of the critical points of R r $ and R r , n , in which the bifurcations are shown. 



Mi 



•*8 i 





■SMI : htJ-. : - 



IS 

***** . Hflatx 



0.4 0.5 0.6 0.7 0.8 0.9 



Asymptotic orbits of critical points of i? ri o (low resolution) . 
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5 



3 



0.4 0.5 0.6 0.7 0.8 0.9 1 0.4 0.5 0.6 0.7 0.8 0.9 1 

Asymptotic orbits of critical points of R r ^ (low resolution). 

In both pictures, there is an obvious complementary nature and there are half-branches 
of bifurcations. These will be explained analytically in the coming propositions. 



Proposition 3.1. The asymmetric 2-cycle of R r ^ bifurcates into a symmetric attracting 
4-cycle at r = 1/y/E. For R r fl, there is a pitch-fork period preserving bifurcation of the 
symmetric 2-cycle into asymmetric attracting ones atr = l/y/2. 

Proof. For the 2-cycle j e l0 , e m j of R r)7T , the derivatives are given by 

1 - 3r 



KM 



1 — r ' 
1 + 3r 



KM 1 + r 



1 - 9r 2 
1 -r 2 

exactly r = 1/ We claim that the new attracting 4-cycle is reflection symmetric. In 



Thus, the 2-cycle undergoes a period-doubling bifurcation when — 1 = — — , which is 



fact, the equation for a reflection symmetric 4-cycle of R r ^ is 

x + ii = l(x) + l o R r n (x) mod tt. 

Let f(x) = x — l(x) — l(x + tt — 2l(x)) + tt. Then we have to solve for f(x) = mod tt. 
Clearly, f(x) — tt is an odd function with /(— tt) = 0, /(0) = tt, and /(tt) = 2tt. If / 
is decreasing at — tt, equivalently, it is so at tt, then / must have a zero modulo 2tt in 
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neighborhoods of ±7r. It is easy to compute that 

/'(-7T) = /'(7T) = 1 - t '(-7T) - .'(o) [i - 2,'H)] 
1 - 5r 2 



1 -r 2 

Thus, /'(— 7i") < if and only if r > l/v5. This shows that for I/a/5 < r, -R rj7r has a 
symmetric 4-cycle. Moreover, it must be attracting when r < 1/ \/5 + £ by continuity 

We have shown in the previous section that an attracting symmetric 2-cycle exists for R r 
with 1/2 < r < l/y/2. It is given by i? rj0 ( c i) — ~ c i mod tt. One may further calculate 
according to the lemma in §2 to obtain that, at r = l/\/2 and x = c±, 

^-^ = 2^-8^0. 
arc's 

Then, by an argument making use of the Inverse Function Theorem, one concludes that 
R r fl has a 2-period preserving pitch fork bifurcation at r = 1/ a/2. □ 



Proposition 3.2. There is a pitch-fork bifurcation on R r 7T where a symmetric orbit of 
period 4 breaks into two asymmetric orbits of period 4- 

Sketch of proof . The key is to consider the zero set of R^^x) = x in the (r, x)-plane. In 
the above, we have shown that for r slightly larger than l/\/5, this can be obtained from 
the solution set of a symmetric 4-cycle, Rl^x) = —x. We then solve for the specific 
values of (r ,x ) such that Rl o 7T (x ) = — x and (i? 2 ^'{xo) = 1. Our calculation involves 
the expression of R r ^ Blaschke product given in 



— 2; 2 (1 — rz) 



z — r 



Then, a symmetric 4-cycle can be solved from z = e lr ^ and 
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After factoring out the obvious factors [z + l)(z — 1), and letting y = cos0, we have the 
polynomial equation 

-1 - 4r 2 + r 4 + 2r(l + 7r 3 )y + 4r 2 (2 - 3r 2 )y 2 - 24rV + 16r 4 y 4 = 0. 

For r > l/y/E, all the four roots of this equation are real. Two of them are always > 1, 
one lies within (—1,0) and the other within (0,1). After taking arccosines, the four 
solutions form a symmetric 4-cycle. Let Xq = Xq{t) be a solution, then we solve for r 
in (i? 2 7r )'(xo) = 1, we obtain a numerical value of r approximately equal to 0.57. We 
further verify that 

- (x ,r )^0. 

Many of the above calculations are lengthy and we indeed make use of computation 
software to help us. Finally, we apply the Inverse Function Theorem to conclude that 
there is a 4-period preserving pitch fork bifurcation. □ 




To end this article, we would like to present two pictures, asymptotic orbits for = 
7r/2, 7r/4 to illustrate the wide variation of the dynamics of this family with respect to Q. 
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Asymptotic orbits of Rr.n/2 an d Rr,-n/A (low resolution). 
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